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Abstract 

We derive the exact general form of the equation of state, in the fu- 
gacity format, for the two-dimensional Coulomb gas. Our results are valid 
in the conducting phase of the Coulomb gas, for temperatures above the 
Kosterlitz-Thouless transition. The derivation of the equation of state 
is based on the knowledge of the general form of the short-distance ex- 
pansion of the correlation functions of the Coulomb gas. We explicitly 
compute the expansion up to order 0(C 6 ) m the activity £. Our results 
are in very good agreement with Monte Carlo simulations at very low 
density. 
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1 Introduction and summary of results 

The system under consideration is a classical two component Coulomb gas com- 
posed of positive and negative particles with charges +1 and — 1. The particles 
live in a two dimensional plane and they are small impenetrable disks of diam- 
eter a. The interaction between two charges q and q' at a distance r from each 
other is 

v(r) = h qq ' ln i r>a (1.1) 
I +oo r < a 

where L is an arbitrary length scale fixing the zero of the potential. This is 
the two dimensional version of the restricted primitive model for electrolytes. 
We shall work using the grand canonical formalism with fugacity A (dimensions 
length^ 2 ) and inverse reduced temperature (coulombic coupling) /?. The arbi- 
trary length scale L can be absorbed in the fugacity by defining the rescaled 
fugacity z = \L^/ 2 which has dimensions lengthy* -4 "" 2 ■ A dimensionless activ- 
ity which will prove useful later can be defined as C = zcr' 4- ' 3 ^ 2 = Acr 2 (L/cr)' 3 / 2 . 
Notice that if L is chosen as L — a, ( — Act 2 does not depend on (3 for fixed 
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fugacity A. In the thermodynamic limit only neutral configurations are relevant. 
Let n + = ri- be the density of positive (negative) particles. The total number 
density is n = n + + ri- = 2n + . 

In the low density limit na 2 — ► 0, there are two values of the coupling of 
special interest. At = 0kt = 4 the system undergoes the Kosterlitz-Thouless 
transition of infinite order [1]. In the high temperature phase < 0kt, the 
system is in a conducting phase with free ions that can screen external charges. 
The correlations have an exponential decay and they satisfy several screening 
sum rules, for instance the Stillingcr-Lovett sum rule [2] . In the low temperature 
phase, for > 0kt, the gas is in a dielectric phase where all charges are bound 
forming dipolar pairs. The perfect screening sum rule is no longer satisfied. 

The other value for of interest is = 2. For < 2 the thermodynamic 
quantities and correlation functions of the system have a finite value in the limit 
of point particles a = 0, while for 2 < < 4 at fixed fugacity z, the density, the 
free energy and internal energy of the system diverge when a — > 0. This is due 
to the collapse of pairs of point particles of opposite sign. On the other hand, 
it is believed [3] that the truncated density correlation functions remain finite 
in the limit er — > when 2 < < 4. 

For < 2 and a — 0, the equation of state for the pressure p of the plasma 
has been known for a long time [4] . A simple scaling argument gives the volume 
dependence of the free energy which leads to the pressure 

/?P=(l-f)n. (1.2) 

On the other hand the temperature dependence of the free energy is highly 
non trivial. Only recently, exact results for the full thermodynamics of the 
two-dimensional Coulomb gas, in the region < 2 and a = 0, have been 
obtained by Samaj and Travenec [5]. These results have been obtained using the 
equivalence between the classical Coulomb gas and the quantum sine-Gordon 
model. In two dimensions this model is integrable, the free energy is known 
in terms of the soliton mass [6], and the relation between the soliton mass 
and the coupling of the sine-Gordon model (i. e. the fugacity of the Coulomb 
gas) in the conformal normalization has been found [7]. This gives the exact 
density-fugacity relationship for the Coulomb gas, which allows to find all the 
thermodynamic quantities of the system [5] . 

In the region 2 < < 4, since the density diverges in the limit no 2 — ► 
for fixed fugacity, it is more appropriate to study the fugacity expansion of the 
pressure, rather than its density expansion. In Ref. [8], Gallavotti and Nicolo 
considered a version of the Coulomb gas with a soft short-distance cutoff. They 
proved that a Mayer series expansion of the pressure in integer powers of the 
fugacity have well defined coefficients up to order 21 for > 0i, where 

l= 4(l-^j 1 = 1,2,3,... (1.3) 

whereas higher order Mayer coefficients diverge. For > 4 all Mayer coefficients 
are finite, while for < 2 all Mayer coefficients diverge. 
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Their findings lead them to conjecture that the plasma undergoes a series 
of intermediate phase transitions at (3 — (3i from the conducting phase at (3 = 
(3i — 2 up to the dielectric phase at (3 = (3^ = /3kt = 4, as opposed to the 
traditional Kosterlitz-Thouless scenario where the conducting-dielectric phase 
transition takes place at (3 = 4. 

Fisher et al. [9] denied this conjecture. They proposed an ansatz for the 
pressure, which, in our notations (£ = zcr( 4 ~^/ 2 ), reads 

oo 

0p = ^(/3)z 4 /< 4 -« [1 + e(f3, ()} + — hiWC ■ (1-4) 

G i=\ 

In this ansatz, they conjectured that b^{(3) and b~2i(P) are analytic for (3 < 4 
and that e(/3, £) is an analytic function of (3 for (3 < 4 and is also analytic in ( 
for C > 0. Furthermore, for no 2 — > 0, at fixed z, e{j3, za^^l 2 ) — » 0. 

While the ansatz (1.4) is fully compatible with Gallavotti and Nicolo findings, 
these later conditions on b^((3), t>2i{(3), and e((3, () imply that the pressure 
exhibits no singularities up to [3 — 4, thus there are no intermediate phase 
transitions. 

Using the exacts results for (3 < 2 and a = [5], Kalinay and Samaj [3] 
devised a method to obtain results for the thermodynamic properties of the 
Coulomb gas in the low density limit ncr 2 <C 1 up to (3 < 3. Their findings 
confirm the form (1.4) of the ansatz proposed by Fisher et al. but the analytic 
structure of the coefficients b^{[3) and b~ii {(3) is different. They have simple 
poles at (3 = (3i but they conjectured that a cancellation occurs. At (3 = (3i, the 
exponent of the fugacity in the nonanalytic part of (3p is integer: 4/ (4 — /3;) = 21. 
Then it turns out that residues of &2/(/3) and b^,((3) at (3 = f3i are opposite, thus 
giving no singularities for the pressure at (3 = (3i, confirming the absence of 
intermediate phase transitions. 

The cancellation of singularities was verified in Ref. [3] at the first threshold 
(3 = (3\ = 2, and conjectured for the other thresholds. The aim of this work 
is to extend further the analysis of Ref. [3]. One important ingredient in the 
analysis of Ref. [3] is the knowledge of the short-distance expansion of the density 
correlations functions of the Coulomb gas, n\_ (r) and n + ' + (r) . The cancellation 
at (3 = 2 was obtained in Ref. [3] using the fact that, at the lowest order when 
r — > 0, n-+_{r) — n+\_(r) <~ z 2 r~ /3 . 

In a recent work [10], we presented the general framework to obtain higher 
order terms of this expansion and explicitly computed the two next order terms 
of the short-distance expansion of the correlation functions. Based on this pre- 
vious analysis [10], we will show that the function e(/3, () in the ansatz (1.4) 
is nonanalytic in (. Actually, we will show that the ansatz (1.4) should be 
generalized to 

00 00 2 

= — E M/?)c 2i +M/*)^ [1 + e i.o(/3, c)]+— Y, T* t 4 -^e n , m (P, C) • 

1=1 m=0 

(1.5) 
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The sum J^* is for n > 2 when m = and for n > for m > 1. For m = 1 
the function e„.i(/3, £) = 0. The functions e„, m and e^o admit an expansion in 
integer powers of £ and In (, 

oo 

e„, m (/9,C) =^g„, m , fc (/3,lnC)C 2fe * + " (1-6) 

oo 

ei l0 (/?,C) = E e iAfe(^ ln OC 2fe+2 (1-7) 

fc=0 

where fc* = fc + 1 for n = 0, 1 and k* = k for n > 2. Notice that the lowest 
order in this expansion is at least ( 2 . Thus, in the limit ner 2 — > 0, fixed z, and 
/3 < 4, these terms are irrelevant in the sense that e„ iTO (/3, zcr^ 4- ' 3 '/ 2 ) — ► 0. 

The explicit calculation of 6^,(/3) and 62 (/?) was done in Ref. [3]. Here we will 
compute explicitly 64 and ei j o,o(/3, In 0- We will show that the cancellation 
mechanism between b^{(3) and &2i(/3) conjectured in Ref. [3] indeed takes place 
for I = 2 at (3 = fo = 3. 

The outline of this paper is the following. In section 2, we will recall some 
basic facts about the exact results [5] for j3 < 2 and a — and about the general 
strategy proposed by Samaj and Kalinay [3] to obtain the thermodynamics of 
the Coulomb gas for (3 > 2 in low density limit. In section 3, we prove that 
the equation of state in the fugacity format has the form proposed in Eq. (1.5). 
In section 4, we compute explicitly the coefficients &4(/3) and ei o,oC$> m C) an d 
verify the cancellation mechanism between b^((3) and b~i((3) at f3 = 02 = 3. In 
section 5, from our analytical results for the equation of state, we compute the 
internal energy and the specific heat of the two-dimensional Coulomb gas and 
we compare our results with ones obtained from Monte Carlo simulations [11]. 



2 Previous results and general strategy 

2.1 The Coulomb gas of point charges for (3 < 2 

For point particles, a — and (3 < 2, the classical Coulomb gas can be mapped 
into the Euclidean quantum sine-Gordon model by carrying out a Hubbard- 
Stratonovich transformation. The grand canonical partition function 5 of the 
Coulomb gas can be written as 



JV(f>exp[-S(z)} 
/Xtyexp[-S(0)] 



with the sine-Gordon action 
where we defined 



S(z) = - J , 



16tt 



0(r)A0(r) + 2zcos(60(r)) 



0/4. 



(2.1) 

(2.2) 
(2.3) 
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Under this mapping, the bulk density and two-body densities of charges q = ±1 
and q' = ±1 are given by [5, 12] 

n q = z(e a ^) (2.4) 

and 

«S(l r - r 'D = z 2 (e ib '" l '^e ibg ' 4 '^) (2.5) 

where the averages are taken with respect to the sine-Gordon action (2.2). 

The complete mapping between the classical Coulomb gas and the sine- 
Gordon model requires [5, 13] to use the conformal normalization [7], where, 
when z — > 0, the free fields are normalized according to (e* b< ^°) e^ lb ^ r ^) z= o = 
r~P. Under this conformal normalization, the expectation value of exponential 
fields is known [14, 15] 



0Q 2 



(e > ^ [WF) ) cMhm (2 ' 6) 

with j(x) = T(x)/r(l — x) where T(x) is the Euler Gamma function, and 



h(Q) 



dt 

T 



sinh 2 (2Q6 2 i) on2 ,,- 2( 



2 sinh(6 2 t) sinh(t) cosh[(l - b 2 )t] 



2Q z b z e 



(2.7) 



This expression is valid for j3\Q\ < 2, otherwise the integral h(Q) diverges, but 
it is possible to do an analytic continuation [16] of this formula for other values 
of (3Q using a reflexion formula satisfied by (e* b( ^) presented in Ref. [15]. For 
Q = ±1, the integral (2.7) can be computed explicitly and 



J ttz \^ ( r^/2) \ tanK/2) 
' J l#/4)j V r (¥)/ ( 4 "^W/4) 1 ' ' 

where £ = /3/(4 - (3). 

Equation (2.4) combined with (2.8) gives the exact density-fugacity rela- 
tionship [5] 

n = n[z,0]=z 4 /( 4 -' 3 ) T !^6 v ,(/3) (2.9) 

with 




W)= (alA „* -nrrr tan«/2). (2.10) 

We introduced the notation n[z, a] to indicate that the density is a function of 
z and the hard core diameter a. Here a — 0, but we shall use that notation 
later on. 

Equation (2.9) together with the thermodynamic relation n = zd((3p)/dz 
leads to 

/3p = b i ,(/3)z^ . (2.11) 
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This is the exact pressurc-fugacity relationship for point-like particles when 
(3 < 2. Notice that both (1.4) and (1.5) are compatible with (2.11) when a = 
and (3 < 2. 

At (3 = 2, the density diverges, as expected, due to the collapse phenomenon. 
As already noticed in Ref. [3], a naive analytic continuation of (2.9) for (3 > 2 
does not give the correct density-fugacity relationship beyond the collapse: the 
function n[z,0] diverges at the thresholds (3 = (3i and can become negative, it 
cannot represent the density in the region 2 < (3 < 4. 



2.2 The Coulomb gas beyond (3 = 2 in the low density 
limit 

The method [3] to study the properties of the Coulomb gas for (3 > 2 and a ^ 
when n 2 a <C 1 is based on the electroneutrality sum rule 



n+ 



= / 

Jm. 2 



n+l_ (r) — n++ (r 



(2) 



dr 



(2.12) 



It is believed that the truncated correlation functions n^ T (r;z,a) are well 



(2,)T 

defined for a = up to (3 < 4. Assuming that the difference n , (r;z,a) — 



rig , (r; z, 0) is negligible for r > a, and using the electroneutrality sum rule (2.12), 
Kalinay and Samaj [3] propose that the density n[z, <r] for the Coulomb gas with 
hard core a is given by 



</</' 



n[z, a] = n[z, 0] — An 



f 

Jo 



n™_(r;z,0)-nf + (r;z,0) 



r dr 



(2.13) 



in the limit no 2 <C 1. 

Using the dominant order term in the small-r expansion of the correlation 



functions, n^l(r) — n++(r) ~ z 



(2) 



-0 



Kalinay and Samaj [3] obtained the first 



correction 



n[z, cr] — n{z, 0] — 47rz z 



2-/3 



(2.14) 



and they showed that it cancels the pole at (3 = 2 from n[z, 0]. 



3 The equation of state in the fugacity format 

3.1 The short distance expansion of correlation functions 
and the operator product expansion 

To proceed further with the program proposed by Kalinay and Samaj [3], we 
need to compute the higher order terms of the short-distance expansion of the 
correlation functions. In a previous work [10], we showed how the operator 
product expansion for the exponential fields in the sine-Gordon model can be 
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used to obtain the short distance expansion of the correlation functions of the 
Coulomb gas. Let us recall and extend some of the results from Ref. [10]. 

The operator product expansion for the exponential fields of the sine-Gordon 
models reads [17] 



n— +00 



^ e ibQ 1 4>(0) e ibQ 2 4>(r)^ _ 



^°Q 2 W<e ib(Ql+Q2+ " )0 > (3.1) 



+C^f Q2 (r)m) 2 (B ( f>) 2 e ib ^ + ^ +n ^) + ■ ■ 

where the dots denote subdominant contributions from higher order descen- 
dant fields l\ t d m *4>]\ 3 8 n i<j) e lb ^ 1+ ^+ n ^, where only the fields with £\ m, = 

rij — m have non vanishing expectation value. The level m = 1 field has 
zero expectation value because it is a total derivative [17]. 

The functions Cq'^ 2 (r) of the operator product expansion have the following 
form [17] 

C n '" 1 (r) = z l™l r m +' 3( 3i ( 32+™/3(Ql+Q2)+2|r l |(l-4)+n 2 /3/2 / -n,m / 2 4-^ ft 2 ] 
Q1Q2 ^ ^ ^ Q1Q2 ^ / V * / 

where each fg^ 2 admit a power series expansion of the form 

00 

fQ?Q,W = T,fZ' m {QiM* k - (3-3) 

fc=0 

The connexion of the operator product expansion with the Coulomb gas is 
clear by noticing that (^KQi+^+n)^ j s closely related to the excess chemical 
potential of an external charge Q\ + Q2 + n introduced in the plasma [18] . Also, 
the coefficients f% ' (Qi, Q2) are expressible in terms (n + 2fc)-fold Coulomb 
type integrals: they are the partition functions of a system with two fixed point 
charges Q\ and Q2 separated by a distance 1 and with n (positive for n > 0, 
negative for n < 0) point particles and k pairs of positive and negative point 
particles approaching the fixed particles. For explicit expressions of some of 
these coefficients see Refs. [17, 10]. 

A few technical details, explained in greater detail in Ref. [10], should be kept 
in mind when using the expansion (3.1) to compute the correlation functions 
of the Coulomb gas. First, the (n + 2fc)-fold Coulomb type integrals in the 
coefficients f^' m (Qi, Q2) are defined for a certain range of values of Qi, Q2 
and (3 since we are dealing with point particles, in order to avoid the collapse. 
Beyond those ranges an analytic continuation should be used. 

Second, for certain values of Q 1 , Q 2 and (3 different terms in the expan- 
sion (3.1) can become of the same order. When this occurs the coefficient of 
each term usually has a pole, but adding all contributions of the same power 
in r and taking the appropriate limit gives a finite result and logarithmic terms 
ln(zr( 4 ~ /3 )/ 2 ) appear. One important case where this happens is in the compu- 
tation of As clearly seen from Eq. (3.1), when Q\ = —Q2 = 1, the terms 
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for n and —n are of the same order in r. One consequence is the appearance of 
hi^zr^ 13 ^ 2 ) terms in the expansion of n+_(r) coming from the contributions of 
terms \n\ > 1 in (3.1). However, for n = 0, the corresponding terms in Eq. (3.1) 
do not have a logarithmic correction. For details see Ref. [10]. This last remark 
is important for the following analysis, and, as we shall see, it is the reason why 
the first term of the equation of state (1.5) is an analytic series in the fugacity. 

Using (2.5) and the operator product expansion (3.1), including the contribu- 
tions of all descendant fields, we find that the general form of the short-distance 
expansion of the correlation functions is 

n= + 0O +OC +00 1 |3(2 + re ) 2 +4m OfcJ-O-l- In I 

-aw = £ EE<i4(--)^^(^) <") 

n= — oo m=0 k=0 

n +-(r) = E LL4vr 2 ) \ zr 2 ) (3 - 5) 

n— — oo m— k— 

The indexes used are m for the order of the descendant field (remember that the 
term for m — 1 is zero), n for the number of particles of sign sgn(n) added and k 
for the number of pair of positive and negative particles added. The coefficients 
n X~m k an d n n m k depend on (3 and eventually on ln(zr( 4 ~ /3 ) / ' 2 ), except those for 
neutral configurations nj|"g k and n+J o k which depend only on (i as explained 
above. 



3.2 The fugacity expansions of the density and the pres- 
sure 

Replacing (3.4) and (3.5) into (2.13) leads to the following form for the density 



-oo +00 -foo 



i[z, v] = n[z, 0]+J2 EE^ c tm,kC^ + <i, fe c 



n— — oo m— k— 



2fc+2+|n 



£2K 



(3.6) 

with n[z,0] given by Eqs. (2.9) and (2.10) and ( = za^ W 2 . Finally, using 
n = zd((3p)/dz, the fugacity expansion of pressure is of the form 

M/3), 4 /( 4 -« 



fip 



(3.7) 



oo +oo +oc 



EE^ 



n— — oo m— A;— 



n.m.k 



2/c+2+|n| 



n — 2,m.k 



2fc+2+|n-2| 



Notice that the terms corresponding to m = and n = are analytic in z. 
These, together with the terms m = and n = 1, reproduce the ansatz (1.4) 
from Fisher ef at. [9]. Writing these terms apart in (3.7) the pressure-fugacity 
relationship can finally be written as announced in the introduction 

oo oo 2 

/?P= ? E M/3)C 2 '+ [1 + ei,oC9, C)] + ^2 E El C^^en.m^, C) ■ 

1=1 m=0 

(3.8) 
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The sum J^* does not contains the terms m = and n = 0, 1 since they are 
explicitly written apart: the sum is for n > 2 when to = and for n > for 
m > 1. 

The functions e n ^ m {0,C) and ei j o(/3, C) admit a power series expansion in 
terms of £ and ln£ given in Eqs. (1.6) and (1.7). They arc at least of order 
( 2 = z 2 a^ 4 ~^ , so they vanish in the limit na 2 — > for fixed z and < 4: they 
are irrelevant. On the other hand, in the analytic part of 0p as a series in ( 
(first sum in (3.8)), the l-th term becomes relevant for > 0i = 4[1 — (1/(2/))], 
i. e. it diverges in the limit na 2 —> 0. 

At p = 0i, the term z 4 /( 4 ~^&,/,(/3) has a pole, as seen from (2.10), but it 
is expected that the coefficient <r~ 2 b2i((3)( 21 , which becomes relevant at that 
point, also has a pole and cancels the divergence from z 4 ^ 4 ^ l3 ^ > b^(0). This was 
checked at = 0\ = 2 for the first term I = 1 in Rcf. [3]. In the next section 
we check that these cancellation also takes place at = 02 = 3. 

4 Explicit calculations 

4.1 The equation of state in the fugacity format 

In Ref. [10] we computed explicitly the short-distance expansion of the correla- 
tion functions (3.4) and (3.5) up to order r 8 ^ 3 ^, that is, we computed the terms 
corresponding to (n,m,k) = (—2,0,0) and (n, m, k) = (—1,0,0) for ra++, and 
(n,m,k) = (0,0,0), (n,m,k) = (0,0,1) and (n,m,k) = (±1,0,0) for n { + ] _. The 
neutral configurations: for n^i, (n, to, k) = (—2, 0, 0) gives the order r 4 ~ 2 ^, and 

for (n,m,k) — (0,0,0) gives the order r^ 13 , and (n,m,k) = (0,0,1) the 

order r 4 ~ 2 ^. The configurations with at most one charge ±1 give the order r 2_/3 
[(n,m,k) = (-1,0,0) for n+| and (n,m,k) = (±1,0,0) for n+1]. Explicitly, 




■V 8 - 3 V 6 - 2/3 ) (4.2) 



r (4.1) 
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with (e^) given by Eq. (2.8), and 



7T/3 2 

4 



-4(l)-4 + 4C 



n. 



+V>(§) + V>(-§) + ^i -§) + </>(! + §) 

T7(l-f) 2 7(/3-l) 

47r2 /, /?x3 / , 3 /?n 

-^5T(l-4) 8 7(-l + T ) 



(4.3) 



(4.4) 
(4.5) 

(4.6) 



where ip(x) = dhiT(x) / dx is the digamma function and C = —tp(l) is the Euler 
constant. The functions I' b (l) — dIb{Q) / OQ\q=\ and J(/3, — /3,/3) are detailed 
in the appendix. 

The use of (2.13) leads to the density 



n[z,a] = z(e a +)i2- 



47rzV-^ 
4-/3 



-+- ~++ 
nJ — n, — 



4-/3 



In 



TTZ \ ' 



4-/3 



- 1 



2-/3 



3-/3 



(n 4 + - 



n+ + ) + OK j( r 



6 ^.10-3/3 ^8-2/3 



Integrating the thermodynamic relation n = zd(j3p)/dz, we find the equation 
of state in the fugacity format 



(4.7) 



1 



0p = bi,(P)zW-l» [l + e(/3,C)] + - 

a 



-2tt 



2-/3 



C 2 + 6 4 (/3)C 4 + 0(C 6 ) 



with 

6 4 (/3) = 

and 

e(/?,C) - 



2(3-/3) 



J(/3,-/3,/3) 



4^ 2 



(2-/3) 



2 7(1 



T) 3 7(-l + 



3/3, 
4 ' 



4^C 2 



(4-/3) 2 (6-/3) 



-^/3 2 



(4.8) 



(4.9) 



(4-/?) 
-^7(l-?) 2 7(/3-l) 



/3, 



- 4 + 4C + V(f ) + ) + ^(1 - |) + V>(1 + |) 



7T/3 2 



2(/3-5) 
6-/3 



+ ln 



7[C A' 
7(0/4) >/ 



+ o(C 2 ) (4.10) 



We confirm that the function e(/3, zcr' 4 ^Z 2 ) — ► when no 1 — > at fixed 
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4.2 Cancellations at j3 — 3 for the relevant terms 

At (3 = 3, the contribution to (3p from n[z, 0] is 

M/?)z 4 /( 4 "« /3 ~ 3 ^ 7 (l/4) 4 z 4 -^ (4.11) 
On the other hand, in the appendix it is shown that 

5^^-^7(1/4)^- (4.12) 

Thus the divergences of each term at /3 = 3 cancel each other yielding a finite 
result for the pressure. The cancellation mechanism conjectured in Ref. [3], 
indeed take place at (3 = 3. 

4.3 Cancellations at j3 — 2 for irrelevant terms 

Actually a more complex mechanism of cancellations of divergences appears 
to take place, at (3 = Pi, when the generically nonanalytic contributions in 
z 4/(4-/3) income analytic z 21 , also with the irrelevant terms (terms that vanish 
when na 2 — > 0). To illustrate this, notice that at P = 2, b^p)^ 4 / 'a 2 has a pole, 
but so does the term b^,{P)z 4 ^ 4 ^ e(p , £) which is also of order £ 4 at (3 = 2. The 
finitcness of the correlation functions in the region up to (3 < 4 (in particular 
at (3 = 2 for the present case) ensures that both divergent contributions cancel 
each other. 

The calculations of the short distance expansion of the correlation functions 
at f3 = 2 has been done in Ref. [10]. Using the results from Ref. [10], in particular 
the results from Appendices A and B of Ref. [10], it is easy to check that the 
contribution of terms of order ( 4 in the pressure are finite: 

7T 3 zV 2 

(3p = 2 -7rz 2 [-l + 2C + 21n(7rzcr)] — { - 3 - 4(C* + ln7r)[l + 2(C* + In 

+4[-l + 21n(zcr) +4(ln7r + C)]ln(zcr)} +o(z 4 cr 2 ) (4. 

We have also written the relevant contribution (nonvanishing when a — > 0) 
which was computed in Ref. [3]. 

If the conjecture that the truncated correlation functions are finite up to 
(3 = 4:, this cancellation mechanism of irrelevant terms should take place at other 
values of (3 where the coefficients in Eq. (3.8) diverge. For instance, in (4.8) the 
product b t f,((3)e((3,() has a pole at (3 = 3. At this value of (3, the nonanalytic 
contribution to the pressure b^,((3)e((3, Qz 4 ^ 4 '^ becomes analytic of order £ 6 . 
The pole at (3 = 3 of b x p((3)e((3 , £) should be canceled with a similar diverging 
term from be(/3)<^ 6 /a 2 , which we have not computed here. 
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5 Comparison with Monte Carlo simulations 



and 



In the canonical format, the excess dimensionless free energy per particle is 
given by 

/(„,/?) = I^P+ In C (5.1) 

where £ should be expressed in terms of the density n by inverting the rela- 
tion (3.6). The excess internal energy per particle u cxc and the excess specific 
heat at constant volume per particle c^? c can be obtained as 

«~ = "-^f- (5.2) 

,« = _, 2 #/M. (53) 

Using the expressions (4.7) and (4.8) for the density and the pressure, accu- 
rate to order 0(C e ), obtained in the previous section, we numerically inverted 
the density-fugacity relationship (4.7) and computed the internal energy and 
specific heat for two low density packing fractions r\ = 7rncr 2 /4, r\ = 5 ■ 10~ 4 
and r\ = 5 • 10~ 3 . In Ref. [11], Monte Carlo simulations of the two-dimensional 
Coulomb gas where performed for these two packing fraction values. 

Our analytical results are compared to the Monte Carlo simulation ones in 
figures 1 and 2. In figure 1, we plot the internal energy u cxc as a function of 
the inverse temperature /?. For rj = 5 • 10 -4 , we show two curves. The dashed 
line corresponds to the results obtained from Eq. (4.8). As it was discussed 
in the preceding section, the term b^((3)e(/3,() from Eq. (4.8) has a pole at 
(3 = 3, which should be canceled with the next order term b 6 ((3)( 6 /a 2 which has 
not been computed. For this reason, the comparison with Monte Carlo results 
can only be done for (3 < 3. However, since this term vanishes when — > 0, 
we decided to compare the Monte Carlo results with those obtained from our 
analytical formulas by omitting this "irrelevant" term. This is shown in figure 1 
with the dot-dash line. Explicitly, this last curve is obtained by approximating 
the pressure by 



(3p 



271 C 2 + h(PK 4 



2-/3 



(5.4) 



instead of using Eq. (4.8). Eq. (5.4) can be used up to the next pole at (3 = 
Ps = 10/3. For very low volume fractions, r\ — 5 • 10~ 4 , the agreement with 
Monte Carlo simulations is very good, both using the correct formula (4.8), for 
(3 < 3, or the "truncated" one (5.4), for 2 < (3 < 10/3. For higher volume 
fractions, r) = 5 • 10~ 3 , the agreement is still very good when using the complete 
formula (4.8). On the other hand, the "truncated" formula Eq. (5.4), docs not 
give a good agreement (curve not shown) as it can expected because the omitted 
term becomes important at high volume fractions. 
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— ■ 11=5 


ID 4 


- n =5 


io- 4 


— T!=5 


ID 3 



1.5 



2.5 



Figure 1: The internal energy u oxc as a function of (3. The lines are our analytical 
results and the bars are the Monte Carlo simulation results. For r\ — 5 • 10~ 4 
the dashed line shows the results obtained with the full expression (4.8), while 
the dash-dot line shows the results obtained from (5.4). 
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I , 1 , 1 , 1 , 1 , 1 , 1 , 1 

1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 

P 

Figure 2: The specific heat c™ c as a function of 0. The lines are our analytical 
results and the bars are the Monte Carlo simulation results. 

In figure 2, we plot the specific heat c^? c as a function of 0. The agreement 
with Monte Carlo simulations is fairly good, even at the relatively high volume 
fraction r] = 5 • 10~ 3 . In any case, the agreement with the simulations is much 
better than the one obtained with only the first order term b 2 (0)( 2 /a 2 in the 
pressure, shown in Fig. 3 of Ref. [3], as expected. 

6 Conclusion 

Using the exact results for the thermodynamics of the two-dimensional Coulomb 
gas of point particles for < 2 [5] , the short-distance expansion of the density 
correlations functions [10], and the program proposed by Kalinay and Samaj [3], 
we have derived the general form of the equation of state in the fugacity for- 
mat (1.5), for the two-dimensional Coulomb gas composed of small core diameter 
a particles and for < 4. We explicitly computed the second corrections due 
to the hard core, up to terms 0(C 6 ) in the activity £, the first corrections (order 
( 2 ) where computed in Ref. [3]. 

The general form of the equation of state (1.5) is compatible with the fact 
that an analytic expansion of the pressure in powers of the fugacity have finite 
Mayer coefficients up to order 21 for > f3i , confirming the findings of Gallavotti 
and Nicolo [8]. However, the explicit calculations we performed show that the 
pressure does not have any singularities at = 0\ = 2 nor at = 02 = 3, 
contrary to the conjecture of a series of intermediate phase transitions proposed 
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by Gallavotti and Nicolo [8] , and supporting the arguments of Fisher et at. [9] 
against this conjecture. But the general form for the equation of state (1.5) 
we found is more complex than the ansatz (1.4) proposed by Fisher et at. [9]. 
Finally we compared our results against Monte Carlo simulations results and 
we found good agreement. 
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A Appendix: Technical details 

The function i£(l) appearing in Eq. (4.3) is defined as i£(l) = dIb{Q) / 8Q\q=\ 
with I b (Q) given by Eq (2.7). Explicitly, 



4e -2t + tsmh(0t) 



sinht cosh[(l - f )i] sinh(/?i/4) 



(A.I) 



This expression converges only for < 2. To use it beyond > 2 it is useful to 
write it as 

1/^(1) =J 1 (/3)+J 2 (/3) (A.2) 

where 



o \ t sinhtcosh[(l - § )t] 

e- 2t cosh(3/3i/4) 
t sinhtcosh[(l - f )t] 



-r{- 



1 2 {0) = 2/| -— + . ; , \[J ' ] dt (A.4) 



The first integral, I\{0), is well defined for < 4. The second one, T 2 (0), is 
defined for /3 < 2, but it can be analytically continued up to < 10/3, by 
writing it as 

. 4 8 2 8 , A „. 

W) - 2-7 + 37^8 + 3-7 + 5^T6 (A ' 5) 



+ 2 



Jo 



e -2t e -t(4+/3)/2 _ e -3t(4-/3) + g-t(4-/3) + g-i(10-3/3) 

~ + ~~ " (l-e- 2 *)(l + e- 2t(1 -^) 



The function J(0, —0, 0) appearing in Eq. (4.4) is a special case of the 
integral 



,' iQ.-HQ.3) = f d '",,J„ ' ( Wi) 



| y |/3Q|l_ a; |/3Q| a ._ y |/3 
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for Q = 1. Using the results from Rcf. [19], in the appendix B of Ref. [10] we 
showed that 



with 



4 



J(f3, -(3,(3) = [s((3/2)J+ + s{P)J+] 2 + [a{fi/2)J+] 

r(i-f) 3 r(2-| 



(A.7) 



r(2 - /3)r(3 - p) 
f)r(2-f; 



( 1-|,2-|.-|:2 - :;- /: I )(A> 



r(i-f)r(2-f)r(i + f) 2 



3^2(2-^,1 



2' 2 
1.00 



T 2 



;2,3;1 .(A.9) 



where we used the notation s(x) = sin(7rx) and where 3-^2(01, 02, 03; 61, 62; 2) 
^ a lli)k(b2)kk\ k zk * s a g encranze d hypergeometric function, and [a)k 
T(a + k)/T(a) is the Pochhammer symbol. 

In the appendix B of Ref. [10] we proved that near (3 = 2, 



J(P,-0,0) = (2tt) 2 

p >2 



1 1 

+ 



09-2)2 p_ 2 



1 



+ 0(J3-2). 



(A.10) 



In order to verify the cancellation mechanism at (3 — 3 presented in Sec. 4.2, 
we need the value of J((3, — /3, (3) at (3 = 3. 

For this purpose, it is convenient to write as 



J+ = 



r(i-f) 3 r(2-§; 



r(4 - (3y 



(2 - /3)(3 - /3) 2 + (1 - |)(2 - f )(-/J/2)(3 - /?) + S(/J) 



where 



<^ = v (l-f)fc(2-f) fc (-/?/2) fc 
lPj ( 4 -/3) fe _ 1 (4-/?) fc _ 2 fc! 



fc=2 



At /3 = 3, we have 5(3) = -7(1/4) 2 /(16tt). Then at /? = 3, 

J+ = %(l/4) 2 . 

On the other hand j£ can be evaluated directly at (3 = 3, using [20] 

3*2(1/2, 5/2, 3/2; 2, 3; 1) = |^7(l/4) 2 . 
We hnd that J+ = - J+ at /3 = 3. This yields 



(A.11) 
(A.12) 

(A.13) 
(A.14) 



J(3,-3,3) = T7 (l/4) 4 . 



(A.15) 



This result, combined with the explicit expression (4.9) for bi{(3), gives the 
behavior (4.12) for b A ((3) near (3 = 3. 
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